Phase transitions and topology in 2 + k XY mean-field models 
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The thermodynamics and topology of mean-field models with 2 + k body interaction terms (gen- 
eralizing XY model) are derived. Focusing on two particular cases (2 + 4 and 2 + 6 body interaction 
terms), a comparison between thermodynamic (phase transition energy, thermodynamically for- 
bidden energy regions) and topological (singularity and curvature of saddle entropy) properties is 
performed. We find that i) a topological change is present at the phase transition energy; however, 
ii) only one topological change occurs, also for those models exhibiting two phase transitions; Hi) 
the order of a phase transition is not completely signaled by the curvature of topological quantities. 



I. INTRODUCTION 



In the recent years different authors suggested a pos- 
sible topological approach to the study of the phase tran- 
sitions. Within this approach it has been suggested that 
any thermodynamic phase transition mirrors a topolog- 
ical change of the potential energy hypersurface [1, 2], 
i.e. a change in the topology of certain submanifolds in 
configuration space (for a recent review see Ref.s [3, 4]). 
More specifically, the energy density e where a thermo- 
dynamic phase transition takes place (e = e c in the mi- 
crocanonical, or e = e(T c ) in the canonical, ensemble) 
has been conjectured to be the same where a topolog- 
ical change of the submanifold M e = {q\V(q) < Ne} 
appears (here q are the generalized coordinates, V{q) the 
potential energy function, N the number of degrees of 
freedom) . This "topological hypothesis" has been subse- 
quently formalized in a theorem which, however, applies 
only to a strict class of systems, described by smooth, 
bounded below, confining, finite range potentials [5-7]. 
This theorem states that a topological change is a nec- 
essary condition for the presence of a phase transition. 
However the study of model systems fulfilling the re- 
quirements of the theorem is a very hard task, due to 
analytical difficulty to solve the thermodynamics and/or 
topology or to numerically calculate topological quanti- 
ties. For this reason, beside few cases [8-10], many works 
have been devoted to the study of tractable model sys- 
tems not-fulfilling the theorem hypotheses [10-16]. For 
such models a variety of results has been obtained, some 
in agreement and some not with the topological hypoth- 
esis (see Table I of Ref . [3] for a summary) . Among the 
others, two particular interesting questions remain to be 
answered. 

The first question concerns the presence of topological 
changes in correspondence of phase transition energy val- 
ues. Is a topological change in those systems undergoing 
a phase transition always present? And more, is the en- 
ergy e T where topological change are observed coincident 
to e c , the thermodynamic phase transition energy? It is 
well established that for certain models (not-fulfilling the 
hypotheses of the theorem) the two energies are not coin- 
cident. Two different mechanisms have been proposed to 



be responsible of this discrepancy: a maximization pro- 
cedure of a smooth function generating a phase transition 
(as opposite mechanism with respect to the topological 
one) [17] or an underlying saddle-dominated dynamics 
for which the relevant topological energy level is not the 
instantaneous potential energy, rather it is the potential 
energy where are located the saddles of V(q) mostly vis- 
ited at the thermodynamic phase transition state point 
("weak" topological hypothesis) [15, 16]. It is worth to 
mention that in all the models analyzed so far there is 
always a topological change (although at an energy not 
coincident with the thermodynamic one) in the presence 
of a phase transition. 

The second open question regards the possibility to 
infer the order of the phase transition from the cur- 
vature properties of topological invariants such as the 
Euler characteristic. In other words, is there a one- 
to-one correspondence between curvatures of thermody- 
namic entropy and some topological invariant quantity? 
In previous studies of toy models (fc-trigonometric model) 
[14], capable of switching between first- and second-order 
phase transition by tuning a control parameter, a posi- 
tive answer to this question was given. Specifically, we 
observed positive curvature of saddle entropy for systems 
undergoing first-order phase transition, while a stan- 
dard negative curvature accompanied the second order 
transitions. The relevant control parameter of the fc- 
trigonometric model is the number fc of interacting bod- 
ies in the Hamiltonian which, in turn, depends on the 
relative phases of these fc-interacting bodies. 

In this paper we study a class of mean-field models with 
different many-body interaction terms, which, therefore, 
can be named as fc+fc' model. Both thermodynamics and 
topology are analytically tractable, so a direct compari- 
son can be made between the two. Two particular model 
systems will be analyzed (2+4 and 2+6), which manifest 
a rich phase diagram. Our main findings are the follow- 
ing, i) There is always a topological change at the same 
energy level for all the model systems, corresponding to 
the paramagnetic energy at which the phase transition 
(first or second order, depending on the model) takes 
place, ii) For one of the two models two phase transitions 
are observed (paramagnetic-magnetic second order tran- 



sition, followed by a magnetic-magnetic first order tran- 
sition), but a topological change is only present in corre- 
spondence of the paramagnetic-magnetic transition. The 
second phase transition seems not to have signature in 
the topology. Hi) A qualitative agreement between sad- 
dle and thermodynamic entropy is obtained (comparing 
positive/negative curvature regions), although the quan- 
titative discrepancy does not allow one to predict the 
presence of a first-order transition from topology for cer- 
tain values of coupling parameters (taking into account 
inherent saddles do not modified the discrepancy, even 
though a better quantitative comparison is obtained). 
Basically, the results of this paper point toward a weaken- 
ing of the link between "thermodynamic" and "topology 
of the potential energy function" in mean-field systems. 
A richer scenario is beginning to appear, and a compre- 
hensive picture is not presently at hand. 

The paper is organized as follow. In Sec. II we will 
introduce the model in its general form, in Sec. Ill wc 
will derive the canonical thermodynamics focusing on two 
particular cases, in Sec. IV the topology will be analyzed, 
calculating stationary point properties. Conclusions will 
be drawn in Sec. V. 



II. THE MODEL 

We consider a class of mean-field Hamiltonians of the 
form 
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FIG. 1: (Color online) H2 + H4 case. Phase diagram in the 
(J4/ J2, PJ2) plane. Second and first order phase transition 
(PT) lines are reported. 



example that will be useful in the following, we explicitly 
write also the k — 3 case 

H& = ~ ^ C ° S ^ + ^ 2 + ^ 3 ~ Vjl ~ Vj2 ~ 

ji ,n ,33 

(5) 

In this paper we will mainly focus on the case of 2- 
terms contributing to Hamiltonian (1) 

H = H 2 +H 2k , (6) 

with k = 2 and k = 3, i.e. H = H 2 +H 4 and H = H 2 +H 6 . 



with 
H 2 k = 



N 2k 



{i},{3} 



-<Ph ' 



(2) 

where the sum is over the sets {i} — ii,...,i k and 
{.]} = Ji, ■ ■ -,3k (h3 = A) and {^} are angu- 

lar variables ifi € [0,27r). We restrict to ferromagnetic 
interactions, i.e. J 2 k > 0. The system described by 
Hamiltonian (1) can be viewed as an ensemble of 2d ro- 
tors interacting through mean-field potentials. For k = 1 
the Hamiltonian (2) reduces to the usual XY mean-field 
Hamiltonian 



H 2 = -^J2 C0S (^ - <Pi) 

i,3 



(3) 



while for k > 1 we have 2fc-body interaction terms. For 
example, for k — 2 



Ha 



J4 \ -> 

31 ,32 



cos ((pi r + ip i2 - ip n - <p h ) (4) 



which has been recently introduced as a model for mode- 
locking laser Hamiltonian [18]. Just to give a further 



III. THERMODYNAMICS 

In this section we derive the canonical thermodynam- 
ics. 

Introducing the complex variable z: 

N 

Y 



(7) 



the Hamiltonian (1) can be written as 



M 



n = -Nj2J2k P 2k 



k=l 



It is now possible to write the partition function 



Z = J d{tp} 



,-pn 



as [191 



Z oc dp e 



(8) 
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FIG. 2: (Color online) H2 + H4 case. Equilibrium magnetiza- 
tion p (upper panel) and energy e (rescaled by J2 + J4, lower 
panel) as a function of inverse temperature j3 (in unit of J2) 
for J4/J2 = 0.1. Second-order phase transition takes place at 
PJ2 = 1. 



where [3 = 1/ksT (ks is the Boltzmann constant, in 
the following we set ks = 1) and the function g(p; (3) is 
explicitly written as 



JY 



A I 



g(p; $) = p 5^(2fc-l) J 2fe p 2fe -ln / 2/3 £ fc J 2fej0 



2fc— 1 



fe=l 



k=l 



(11) 

where / Q (a:) is the modified Bessel function of order 0. 

Performing the thermodynamic limit N — ► 00, the 
saddle-point solution dominates the integral in (10). The 
saddle point equation is written as: 



P 



I (2f3j: k kJ 2k p^) 



(12) 



This equation can have many solutions for p, that with 
the lowest free energy is the stable one 



PfeqiP) = min#(p;/3) . 
p 



(13) 



Thermodynamic properties are obtained from the (3- 
dependence of g. We note that the "paramagnetic" 
(p = 0) solutions always exist and is stable for small 
(3 value. On increasing f3 solutions with p ^ becomes 
possible and, eventually, stable. 

We note that for J 2 = only first order phase tran- 
sitions are possible. Indeed, in order to have a second 
order transition the curvature of the free energy, or of 
the function g, has to change around the paramagnetic 
solution p = 0, from positive to negative. For J 2 > 0, 
close to p = 0, we can expand Eq. (12) obtaining: 



g ~/3J 2 (1-/372)^ 



(14) 



then a second-order phase transition takes place at (3 J 2 = 
1 (if not prevented by a first-order phase transition at 



FIG. 3: (Color online) H2 + H± case. Equilibrium (full lines) 
and metastable (dashed lines) magnetization p (upper panel) 
and energy e (rescaled by J2 + J4, lower panel) as a function 
of inverse temperature /3 (in unit of J2) for J4/J2 = 1.0. A 
first-order phase transition takes place at /3J2 — 0.77. ei 
is the energy corresponding to the appearance of metastable 
free-energy states. 



higher temperature). For J 2 = instead, denoting with 
J p the first non-zero term (p > 4), the expansion of Eq. 
(12) reads: 



g ~ f3J p p p , 



(15) 



the curvature is always positive and then second-order 
transition cannot take place. 

For a single term Hamiltonian TL = H2k, the sys- 
tem undergoes second- (k = 1) or first- (k > 2) order 
phase transition. We do not analyze these cases here, as 
the thermodynamic/topology relationship falls into the 
same class of similar model systems - XY [10] and k- 
trigonometric [14] already discussed in the literature. We 
focus our attention in many-terms Hamiltonian, taking 
into account two interesting cases. 



A. H2 + H4 case 

Considering the Hamiltonian TL = H2 + i?4, similarly 
to the corresponding Ising case [20], the plane spanned 
by the two coupling parameters is split in two regions 
(see Fig. 1), which, as can be seen from Eq. 12, are only 
determined by the ratio J4/ J 2 . Specifically, 

i) For J4/J2 < 1/4 a second order phase transition 
takes place at /3J 2 = 1. As an example, in Fig. 2 the 
equilibrium magnetization p (upper panel) and energy 
per particle e = — J 2 p 2 — J^p (normalized by J 2 + J4, 
lower panel) are reported as a function of the (scaled) 
inverse temperature /3J 2 for the specific case J4/ J 2 = 0.1. 

ii) For J4/J2 > 1/4 the transition becomes first or- 
der, with a jump in both magnetization and energy - full 



4 



1.0 



CO. 



0.5 



II" order PT line 

/ 



I" order PT line 

i 



0.5 



1.0 
J 6 /J 2 



1.5 



2.0 



FIG. 4: (Color online) H2 + He case. Phase diagram in the 
(Je/ J2, 0J2) plane. Second and first order phase transition 
(PT) lines are reported. 



lines in Fig. 3 where the same quantities as in Fig. 2 are 
reported for the specific case J4/J2 = 1-0. The dashed 
lines in Fig. 3 represent metastable states: metastable 
minimum of free energy for e < e\ and local maximum 
for e > e%. 

For the particular case J4/J2 = 1/4 a tricritical point is 
present at PJ2 = 1. 



B. H2 + He case 

The case H = H2 + H§ has a richer phenomenology 
(see Fig. 4). 

i) For Jq/J2 < 0.46 the system, similarly to the case 
H2 + Hi, undergoes a second order phase transition at 
(3 J2 = 1. Fig. 5 reports, as an example, the order parame- 
ter and the energy as a function of the scaled temperature 
for the specific case J§/ J2 =0.1. 

ii) For 0.46 < Jq/J2 < 0.66 a first order phase transi- 
tion occurs after (on increasing j3, lowering T) the second 
order one (Fig. 6 reports the specific case Jq/ J2 = 0.6). 

Hi) For Jq/J2 > 0.66 only the first order transition 
survives (Fig. 7 reports the specific case J&/ J2 — 1-0). 

The dashed lines in Fig.s 6 and 7 represent metastable 
states: metastable minima for e < e± and e > e2 and 
local maximum for e\ < e < 62. 



IV. TOPOLOGY 

In this Section we will study the topology of the two 
models introduced in the previous Section. 

We will follow the same line of calculation performed 
on similar models in previous works [14, 18, 21]. A quan- 
tity directly related to the Euler characteristic of the 
manifold M e = {ip\H((p) < Ne}, is the configurational 
entropy of saddles [14]: 




FIG. 5: (Color online) H2 + He case. Equilibrium magneti- 
zation p (upper panel) and energy e (normalized by J2 + Je, 
lower panel) as a function of inverse temperature /3 (in unit 
of J2) for Je/J2 = 0.1. Second-order phase transition takes 
place at PJ2 = 1. 



where n(e) is the fractional saddle order (i.e. the fraction 
of negative curvatures at the saddle points that are found 
in the potential energy hyper-surface when V(q) = Ne). 
Indeed, we can make the following arguments. 
A stationary point is defined by 



OH 



M 



2kJ 2k p 2 "- 1 ninfo - V) = Vj (17) 



fe=i 



The solutions with p > (those with p = are located 
at e = energy) are obtained from sm(<^j — ip) = (for 
all j), then 



tpj = bP + mjTr] mod 27r 



(18) 



where rrij = {0, 1}. Substituting this solution in Eq. (7) 
we obtain 



P : 



where 



1 -In 



N ^ 



(19) 



(20) 



is the fractional saddle order (as it will be clear soon) 
and we have used the identity (— l) m i = 1 — 2mj. 
We conclude that there are no stationary points with 
n > 1/2 (p is positive defined). The order of a station- 
ary point is defined by its downward curvatures, i.e. by 
the number of negative eigenvalues of the Hessian matrix 
Hij = d 2 H /difiidipj. It is possible to show that in the 
thermodynamic limit the Hessian becomes diagonal 



cr(e) = —n(e) lnn(e) — (1 — n(e)) ln(l — n(e)) , (16) 



H, 



(21) 
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FIG. 6: (Color online) Hi + He case. Equilibrium (full lines) 
and metastable (dashed lines) magnetization p (upper panel) 
and energy e (normalized by J2 + Je, lower panel) as a function 
of inverse temperature /3 (in unit of J2) for the specific case 
Je/ J2 = 0.6. A second-order phase transition takes place at 
PJ2 = 1, followed by a first-order one at PJ2 — 1.02. ei 
and ei are the energies corresponding to the appearance of 
metastable free-energy states. 



FIG. 7: (Color online) Hi + He case. Equilibrium (full lines) 
and metastable (dashed lines) magnetization p (upper panel) 
and energy e (normalized by Ji + Je, lower panel) as a function 
of inverse temperature /3 (in unit of J2) for the specific case 
J&/J2 = 1.0. A first-order phase transition takes place at 
PJ2 — 0.88. ei and ei are the energies corresponding to the 
appearance of metastable free-energy states. 



where 



M 



Xj = cos(<£j - ip) 2fc J 2 fe p 2k 1 



fc=i 



A I 



(22) 



(-i) mj J2 2kJ - 



2k P 



„2fe-l 



fc=l 



Therefore, the saddle order is given by the number of 
rrij = 1 at the considered saddle point, and the fractional 
saddle order n is then given by Eq. (20). Moreover, the 
number of saddles with a given n is given by the binomial 
coefficient, then 



1 f N 
a(e) = lim — In „ , , . 
v ' n^oo N \Nn(e)J 



(23) 



from which the Eq. (16) follows straightforward. The 
latter can be written in the form 

.(e) = -^-M In ±=M - i±^M ln l±M 
w 2 2 2 2 

(24) 

where p(e) is obtained from the thermodynamics, i.e. 
from the solution of the equation: 



e = -^J 2fe p 2k {e) 



(25) 



k=l 



We note that the quantity a is singular at e — 0, due to 
the fact that e > is a forbidden energy region, so a has 
a discontinuity, jumping from a finite value to zero. 



As discussed in the introduction, a central quantity in 
the comparison between thermodynamic and topology is 
the curvature of cr(e). Specifically, we are interested in 
finding the energy values where there is changes in curva- 
ture of cr(e), i.e. the energies where the second derivative 
of er(e) vanishes. After some algebra we get 



d 2 



1 (dp 



T7 = "o X S IP( 



(26) 



where the function s, which depends on e through p(e), 
is: 



s(p) 



1 



1 



1-p l+p 
From Eq. 26 we have 



J2 k 2fc(2fc - 1) J 2k p 2k - 2 ^ 1-p 



s I 



de 2 < 



,2/s-l 



1 + P 

(27) 



(28) 



that is upward, null, downward curvature respectively. 
Studying the positivity of s allows us to determine the 
curvature of a. We now specialize the calculations to the 
two previous cases. 



A. H2 + H4, case 

In Fig. 8, upper panel, the quantity a(e) is plotted as a 
function of energy e (normalized by J2 + J4) for three se- 
lected case J4/J2 = 0.1,0.4,1.5. Full lines correspond to 
negative curvature, while dashed lines to positive (sym- 
bols represent turning points). In the J4/J2 — 0.1 case 
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FIG. 8: (Color online) H2 + H4 case, (a) Saddle-entropy 
a as a function of energy e (normalized by J2 + J4) for 
J4/J2 = 0.1, 0.4, 1.5. Full (dashed) lines correspond to nega- 
tive (positive) curvature. Symbols mark turning points, (b) 
Couplings-energy plane: light-grey region corresponds to a 
positive curvature, dark-grey region to entropy S positive cur- 
vature. Dashed line is the inherent saddle counterpart of the 
border of dark region. 



FIG. 9: (Color online) H2 + H& case, (a) Saddle-entropy a as 
a function of energy e (normalized by J2 + Je) for the specific 
cases J&/J2 = 0.1,0.4,1.5. Full (dashed) lines correspond to 
negative (positive) curvature. Symbols mark turning points, 
(b) Couplings-energy plane: light-grey region corresponds to 
a positive curvature, dark-grey region to entropy S positive 
curvature. Dashed line is the inherent saddle counterpart of 
the border of the dark-grey region. 



(second order phase transition located at e = 0) the cur- 
vature is always negative, while in the other two both 
regions are present. The quantity cr(e) is singular (dis- 
continuous) at e = 0, corresponding to the thermody- 
namic transition energy. In Fig. 8, lower panel, the plane 
( Ji/ J2, e) is drawn (the energy is normalized by J2 + J4). 
The light-grey region corresponds to positive curvature 
of cr(e), its border being the null curvature line - the cor- 
respondence with the turning points of tr(e) in the upper 
panel of Fig. 8 is evidenced with the thin dashed lines for 
J4/J2 — 0.4 and 1.5. The border of the light-grey region 
does not intersect the value J4/J2 = 0.1, and, therefore, 
no turning point exists for such a J4/ J2 value. 

The dark-grey region (that fully lies inside the light- 
grey one) represents the thermodynamically forbidden re- 
gion e > ei (see Fig. 3). It corresponds to positive curva- 
ture of thermodynamic entropy 5(e). If the topological 
hypothesis was correct, the two, light- and dark-grey, re- 
gions would coincide. This is not the case, suggesting a 
not one-to-one correspondence between thermodynamic 
and saddle entropy. 

In the study of other models [15, 16], when the en- 
ergy of the topological change were not coincident with 
the energy of the phase transition, it was found that the 
"weak" topological hypothesis was correct. Indeed, it has 
been shown that the correspondence between topological 
change and phase transition energies were obtained con- 
sidering inherent saddle properties: the energy of topol- 
ogy transition has been found to correspond to the in- 
herent saddle energy. The latter quantity was obtained 
minimizing the quantity W = \W\ 2 [16, 21, 22]. 



Dashed line in Fig. 8b is the inherent saddle line ob- 
tained from the border of the dark-grey region (zero cur- 
vature thermodynamic entropy) applying the minimiza- 
tion procedure of W. The latter has been performed 
solving steepest descent equations of the form <p = — VW 
starting from equilibrium initial configurations and tak- 
ing at the end the infinite time limit. Following similar 
calculations of Ref. [21], we obtain the saddle energy 
e s = e(poo) where poo = L (—f3de/dp)/Io(-(3de/dp) and 
La is the modified Struve function of order 0: Lo(a) — 

2ir~ 1 J^ 2 dip smh(a cos ip). A better correspondence is 
obtained between the light-grey border and the dashed 
line, however the two regions are not yet coincident. This 
indicates that the "weak" topological hypothesis could 
only be considered a good approximation, but not a 
quantitative prescription for the location of the phase 
transitions. A comment is in order. Although obtained 
from canonical ensemble, the saddle energy map (from in- 
stantaneous to saddle energy) is expected to be ensemble- 
independent, when extrapolated to thermodinamically 
forbidden energy regions. "Ensemble inequivalence" phe- 
nomena [23] are then supposed to not affect the results. 

Summarizing: the singular behavior of a at e = sig- 
nals the presence of a thermodynamic transition, the cur- 
vature of a is "quasi" -related to the presence of thermo- 
dynamically forbidden region and then to the appearance 
of a first-order transition. The "quasi" -relation is due to 
the fact that there are regions where the curvature of a 
is positive and the curvature of thermodynamic entropy 
is negative. 
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B. H2 + He, case 

In Fig. 9 the same as in the previous case is reported for 
Ti = H 2 + H e . In the upper panel the energy dependence 
of a is reported for the specific cases J 6 / J 2 = 0.1, 0.4, 1.5. 
As before, dashed lines correspond to positive curvature 
regions. In Fig. 9b the plane (Je/J2,e) is drawn: light- 
grey region is the cr-positive curvature, dark-grey one is 
the S'-positive curvature. Dashed line is the inherent sad- 
dle counterpart of dark region border. Again, beside an 
overall qualitative behavior, we do not find a quantita- 
tive correspondence between the two regions that, ac- 
cording to the weak topological hypothesis, should co- 
incide. We note that there are values of parameters 
(0.2 < Jq/ J 2 < 0.5) where the curvature of a is upward 
(for certain energy values) and a second order transition 
takes place (the curvature of S is downward). 

V. CONCLUSIONS 

By analytically studying the thermodynamics and the 
topology of mean-field models obtained by the sum of 
two different interaction terms in the Hamiltonian (1) - 
2 + 4 and 2 + 6 body terms - we are able to compare 
thermodynamic and topological quantities and test the 
validity of both the "topological hypothesis" [1] and the 
"weak topological hypothesis" [16]. The models have a 
rich phase space structure. The H 2 + H4 model performs 
second or first order phase transition depending on the 
coupling parameters values (a tricritical point joins the 
two) . The H2 + Hq model has a further possible behav- 
ior: it can undergo a double phase transition, a second 
order followed by a first order one. Topological invariant 
(saddle entropy) has a change in correspondence of the 



zero energy (the paramagnetic energy), so confirming in 
that case the equivalence between topological and ther- 
modynamic transition points. However, only one topo- 
logical change occurs, also for the H 2 + H 6 model ex- 
hibiting two phase transition points. Indeed, in this case 
the first order phase transition does not have a topolog- 
ical counterpart, being the topological saddle entropy a 
smooth function at the corresponding phase transition 
energy values. It is worth noting that this is a quite un- 
expected result, all the models analyzed so far presenting 
a topological change at some energy value (even though 
not coincident with the phase transition one for some 
model system). Future studies should establish if this is 
a mean-field "pathology" or has a deeper origin. More- 
over, the curvature of the saddle entropy (as a function of 
energy), seems not to be strictly related to the presence 
of first order phase transition (as previously observed in 
a different model [14]): there are regions in the param- 
eters values for which positive curvature corresponds to 
second order phase transitions. In other words, the cur- 
vature properties of saddle entropy do not coincide with 
those of the thermodynamic entropy. Taking into ac- 
count the possibility that the relevant energy levels are 
given by underlying saddles, a better agreement between 
curvature regions has been found, although there are no 
quantitative coincidence. 

In conclusion, for the analyzed mean-field models a 
topology change is present at the same energy level at 
which a phase transition takes place, in agreement with 
the "topological hypothesis". However, the information 
encoded in the topology seems not to be sufficient to 
predict all the possible thermodynamic behaviors of the 
system, like the presence of two phase transition points 
or the phase transition order. 
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